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ABSTRACT 


This work was undertaken to determine the pressure drops in the 
oxidizer flow passages of a compact dual-combustion liquid rocket 
engine by means of fluid dynamic model tests, using air instead of 
liquids, to establish design data. Investigated also were flow dis- 
tributions and flow instabilities by measurements and flow visuali- 
zation methods. Changes are indicated to improve the original design. 

The air tests were carried out over a range of Reynolds numbers 
to establish the effect of compressibility on the pressure drop loss 
coefficients. Water tests were undertaken to check the validity of 
the results and of the proposed calculating method. 

The tests were conducted at the Turbo-Propulsion Laboratory of 


the Naval Postgraduate School, Monterey, California. 


Thesis by Robert Joseph Kelly entitled: 


of the Oxidizer Flow Passages in a Dual- 


Rocket Engine" 
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SECTION 1 


INTRODUCTION 


The science of astronautics has progressed at an exponential rate 
in the last decade. The development of liquid rocket engines capable 
of producing large thrusts at high specific impulse has been responsi- 
ble for this progress. These engines require very large flow rates of 
fuel and oxidizer. To keep their size and weight at a minimum the com- 
bustion pressures are high, limited only by structural and turbo-pump 
design considerations. The flow passages inside the units are intri- 
cate and costly to manufacture. Pressure losses in these passages are 
high and must be known in order to establish pump and turbine design 
point data. Accurate prediction of these losses by analytical means is 
not possible except for flow passages with simple geometries. There- 
fore, it is advisable to conduct model tests with fluids at low pres- 
sure, preferably air, to determine pressure loss coefficients experi- 
mentally which can then be used to predict the pressure losses in the 
actual engine. 

The ARES (Advanced Rocket Engine Storable) propulsion unit is a 
compact dual-combustion liquid rocket engine designed by the Aerojet- 
General Corporation, Sacramento, California. Professor M. H. Vavra of 
the Department of Aeronautics, Naval Postgraduate School, Monterey, 
California, has goat ted full-scale fluid dynamics models of this unit 
and has conducted preliminary tests of the oxidizer flow passages with 
air. The models have been made available to the writer by Aerojet- 
General and were used to carry out more extensive and more detailed 


testing with improved instrumentation. 
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The purpose of this thesis is threefold. First, the air test 
results are used to predict the pressure drops in the flow passages, 
flow distribution within the unit, and velocity distribution at the 
entrance to the combustion chamber for a liquid flow. Second, the 
effects of compressibility and Reynolds number on the test results are 
investigated, and the validity of the testing method is established. 
Third, changes in design are recommended based on the test results and 
flow visualizations. 

The author wishes to acknowledge the guidance provided by Profes- 
sor Vavra and extends sincere thanks to Mr. R. W. Savage for his assis- 
tance in performing experiments. Acknowledgment is due also to Mr. R. 
Beichel, Aerojet-General Corporation, for the loan of the models and 
his endeavors to apply model testing techniques to the design of ad- 


vanced rocket engines. 
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SECTION 2 


DESCRIPTION OF ARES UNIT AND TEST MODEL 


A three-dimensional schematic representation of the oxidizer 
passages of the ARES propulsion unit is shown in Fig. 1. The oxidizer 
is composed of a mixture of nitrogen tetroxide and hydrazine and it is 
used for regenerative cooling of the combustion chamber and the rocket 
nozzle before being mixed with the fuel. Oxidizer enters the outer 
annulus through a centrifugal pump. This annulus contains turning 
vanes, not shown in Fig. 1, to distribute the flow uniformly to the 
tubes that form the outer skin of the nozzle. Each tube is an inte- 
gral unit consisting of a down-flow and an up-flow passage., The en- 
trance and exit holes are in the same plane perpendicular to the nozzle 
axis; however, they are located at different radii in this plane. The 
tubes begin and end with an inner diameter of 0.3 inches. At the noz- 
zle throat they form an ellipse with major and minor axes of 0.31 in. 
and 0.19 in., respectively. This elliptical cross section gradually 
increases until at the bottom of the tubes the major and minor axes 
have values of 0.76 in. and 0.70 in. A scale model of one of these 
tubes is shown in Fig. 2. From the discharge of the tubes the oxidizer 
enters a second annular chamber concentric with that into which the 
pump discharges. In this annular passage the flow is turned to enter 
a series of 60 holes drilled radially in the inner shell of the second 
annulus. The oxidizer is then discharged into the primary combustion 
chamber through 180 injectors. A flow passage, not shown in Fig. l, 
directs the combustion gases to a single-stage axial turbine which 


drives the pump previously mentioned. Details of this flow passage are 
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given by eee Tests of two turbine designs have been performed, the 
results of which are reported by serena and Messegee. 

The portion of the unit from the inlet of the outer annulus to the 
bottom of the rocket nozzle is called the oxidizer discharge. The por- 
tion including the remainder of the nozzle tubes to the injectors is 
called the oxidizer return. Secondary combustion takes place after the 
turbine discharge and ahead of the nozzle throat. 

A scale model of these flow passages has several modifications as 
far as shape and layout of the actual unit is concerned, but mostly re- 
tains flow passage dimensions. The model, shown assembled in Fig. 3, 
is constructed of plastic and aluminum and weighs 136 pounds. A sche- 
matic of the flow passages is shown in Fig. 4. The discharge and re- 
turn annuli are separated by 7 parts of varying flow areas and shapes. 
Flow enters the model through simulated pump vanes which impart a whirl 
component. The whirl is removed by turning vanes and the flow is dis- 
tributed to part 1, shown separately in Fig. 5. Part 1 has 32 equally 
spaced groups of 3 holes with a total flow area of 4.712 sq. in. Parts 
2 and 3 are shown respectively in Figs. 6 and 7. Each of these parts 
has 52 equally spaced single holes with a total flow area of 2.552 sq. 
in. Parts 4, 5, and 6 together with 2 orifice plates simulate the noz- 
zle tubes. An exploded view of this assembly is shown in Fig. 8. The 
first orifice plate, having holes of diameter d = 0.217 in. and dia- 
meter ratio d/D = 6 = 0.723 simulates the downhtftow ‘ofi.itha nozzle 
tubes. The second plate with d = 0.216 in. and B = 0.720 simulates 
the up-flow of the nozzle tubes. Part 7, shown in Fig. 9, has 32 
equally spaced triple holes with a total flow area of 5.574 sq. in. 


From part 7 the flow passes into the second annulus shown in Fig. 10, 
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where it is turned to enter the injectors. Each of the 60 injector 
tubes has 3 holes labeled A, B, C, respectively, from the outer radius 
inward. The tubes are labeled clockwise from 1 to 60 with the first 
tube just past the nine o'clock position when the injectors are viewed 
from the combustion chamber (see Fig. 36). 

The model can be divided into sub-assemblies for testing. Tests 
on the oxidizer discharge sub-assembly including the first annulus 
through part 5, excluding the second orifice plate, were made by vee. 
Preliminary tests of the return sub-assembly, part 6 through the in- 
jectors including the second orifice plate, were also made by Vavra. 
This paper presents more complete tests of the return sub-assembly, in 


addition to tests of the fully assembled unit. 
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SECTION 3 


TEST INSTALLATION AND INSTRUMENTATION 


Available as a source of air is a 200 H. P. centrifugal compres- 
sor which delivers about 4 pounds mass of air per second at a pressure 
ratio of about 2. Part of this flow passes through a cooler; the rest 
discharges through a butterfly valve to atmosphere. Figure 13 is a 
schematic of the flow measuring installation. The cooled air enters 
plenum 1 then passes through a 4 in. pipe discharging into plenum 2 
shown in Fig. 1l. Flow measurement is accomplished by either a sharp 
edge orifice or a flow nozzle installed in the 4 in. pipe (Fig. 12). 
The design of the orifice and the flow nozzle is in accordance with VDI 
aeanderndene and the installation satisfies or exceeds the specified re- 
quirements of pipe lengths upstream and downstream of it. In addition 
a flow straightener is installed at the discharge of plenum 1. The 
Orifice has a diameter d = 3.024 in. and a diameter ratio 6 = 0.7429. 
The nozzle has d # 3.022 in. and 8 = 0.734. Initial tests conducted 
by Vavra used the sharp edge orifice to measure flow rate. Prior to all 
tests performed by the writer the flow nozzle was installed. 

Appendix A contains the derivation of the flow measuring equation 
for the nozzle. Similar information on the orifice is given by Kelly. 
Plastic hoses connect pressure taps with the measuring instruments. 


The pressure difference Ps - P, is measured by a Meriam Vernier Mano- 


d 
meter filled with water. The gage pressure ahead of the nozzle is 
measured by a similar manometer filled with mercury. The measuring 
accuracy of these manometers is about + 0/002. in. The .temperature t 


: Oo 
is measured with an accuracy of + 1 F. by an iron-constantan thermo- 


couple located 48.5 in. upstream of the nozzle. 
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Within plenum 2 there is arranged a bend and a 10 deg. conical dif- 
fuser with a series of screens at the discharge to minimize flow insta- 
bilities. A second thermocouple is situated within plenum 2. Plenum 
pressure is read on either a mercury filled manometer or a water filled 
manometer bank. Three such manometer banks were used. Bank 1 is a 
250 cm. 40-tube bank which can be read to an accuracy of + 0.1 cm. 

Bank 2 is a 175 cm. 60-tube bank with a reading accuracy of + 0.2 cm. 
Bank 3 is a 92 in. 20-tube board that can be read with an accuracy of 
+ 0.1 in. 

The model is attached to plenum 2 by means of a steel plate visi- 
ble in Fig. 11. Mounted on this plate is a wooden nozzle of circular 
cross section to produce uniform entrance conditions. Situated at the 
throat of the nozzle is a Kiel probe measuring total pressure aoe 
For tests with the entire model assembly the conduit shown in Fig. 14 
is used to introduce air from opposite directions into the first annu- 
lus. Static pressure taps measuring Po are located 1.5 in. Blead of 
the model inlet. In the same plane are located probes to measure P. ‘ 
One such probe is pictured in Fig. 15 and shown schematically in : 
Fig. 16. The design allows rotation of the probe into the flow, en- 
Suring an accurate total pressure reading. 

Figure 4 shows the location of the measuring planes within the 
model. Pressure taps measure static or total pressure depending on 
their location. In each of the planes Z, B, C and E are arranged 
8 equally spaced static taps. Planes Y and G each have 8 equally 
Spaced total taps. Planes A and F have 4 static and 4 total taps, 


spaced alternately. Plane H has 12 equally spaced static taps. Plane 


X contains no taps, but the pressure can be accurately predicted at 
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this location. Total and static pressures are measured at selected 
injector nozzles in plane I by the probes shown in Fig. 17. These 
probes, having a 15-deg. chamfer at the inlet, were designed to mea- 
sure total pressure of flows with yaw angles up to 20 deg. From 
Vavra's results it was suspected that the flow did not leave each in- 
jector axially. Therefore the new probe assembly was designed to per- 
mit quick adjustment to reach the condition depicted in Fig. 18. Total 
and static pressures at seven locations within the second annulus were 
obtained by inserting a removable probe which could be oriented into 
the direction of flow. All pressure taps were connected to the mano- 
meter bank by plastic tubing and measured against atmospheric or ple- 
num pressure depending on the test configuration. 

Testing of the return sub-assembly required a modification of the 
installation. An adapter fitted with screens, shown in Fig. 20, was 
mounted on plenum 2. Attached between the adapter and the sub- 
assembly was a ring of 52 inlet nozzles to introduce the flow uniformly 
to the orifice ring. Plane D, only used with this configuration, has 4 
pressure taps. Because flow is essentially stagnant within the adapter, 


it is assumed that PD = P. . A view of this test installation, without 
D 


the pressure taps connected, is shown in Fig. 21. 
The inlet pressure for both configurations was remarkably steady 
and, except where noted, the pressures within the model did not fluc- 


tuate more than + 0.5 cm. of water even at the highest flow rates. 
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SECTION 4 


THEORETICAL CONSIDERATIONS 


If a steady incompressible, adiabatic, viscous flow exists in a 
pipe, Bernouillis equation between two points on a streamline can be 


written as 


2 2 
oa Vo 
(144)P, + € a + ez, = (144)P, + Ss +ez,+h | (1) 
where: 
P = static pressure Guine 
‘ 3 
@ = density (1b /ft”) 


V * velocity (ft/sec) 


a2 7 (Eiieee> 


09 
| 


Z 


x 


height above reference (ft) 


head loss Ciny te) 


Since the first three terms on each side of equation (1) represent 


the total pressure, the head loss is equal to the total pressure loss 


AP in the pipe. 
1-2 
Experience shows that this loss can be expressed by the function 


oP =fCe .K,V,4, 4, e) (2) 
1-2 
where: M ™ viscosity (Cibwee/ fe) 
or = characteristic length (ft) 
d = characteristic diameter (ft) 
e = roughness parameter based on average variation of inside 
radius (ft) 
Buckingham's 7 theorem’ states that the number of independent 


dimensionless groups used to describe a phenomenon with n vari- 


ables is equal to n-r where r is the number of basic dimensions. For 


Zo 


the basic dimensions of force F, length L, and time T, a dimensional 


analysis of equation (2) gives 


Ap 

1-2 S rc ane 
a = f 3 wena d (3) 
e Vv /g BH 


Equation (3) has four independent parameters 


Ap 
1-2 
= ss Euler number 
ev /g 
e vd 
“EK = Reynolds number = Re 
7" relative roughness 
 - 
A geometrical factor 


Usually equation (3) is expressed by 


2 
a, 7 54 a ) 
where S is the friction loss coefficient that depends on the relative 
roughness and Re. In pipe flow with Re higher than about 10° the 
coefficient S remains essentially constant. In addition, the so- 
called momentum loss occurs if pipes have either sudden contractions 
or expansions. It is primarily due to eddying turbulence and has been 


found experimentally to satisfy the relation 


2 
AP V 
tio = ni oe (5) 


where dr is the momentum loss coefficient which depends primarily on 


flow geometry and remains practically independent of Re. 
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In a flow channel where frictional and momentum losses occur, the 


two loss coefficients may be combined to give 


2 

te, . 1 Ss 4 (6) 
where W is a total loss coefficient which depends on Re. If Re is 
sufficiently high, the coefficient v is a function of geometry 
alone. 

The analysis can be extended to systems of flow passages by con- 
sidering an equivalent pipe for which equation (6) is valid. If such 
a system is modeled to scale, the factors YW and £ are contained in 
one loss coefficient 5 . For air as a testing fluid, care must be 
taken to maintain low Mach numbers M in the model as equation (6) holds 
for incompressible flows only. 


For air at low M and with the velocity V, at point 1 as a reference 


1 


vi Pow (7) 


el 


where: w = flow rate (1b _/sec) 





A, ™ area «Ge» 


Equation (6) may then be written as 


2 
AP 4 W 
ty. 7% Sree 7 (8) 
7 e Ay 


Using the equation of state for a perfect gas 


p 
_ 1 (144) 
C "Rr >) 
t 


i 


Z3 


£t-1b 


where R = 53.35 ; 
1b. R 


T, = total temperature 
1 


equation (8) becomes 


t 
1-2 = l 5 
P 2 
ty 2A, 1-2 


(°R) 


(10) 


The so-called equivalent flow rate referred to total conditions at 


1 is defined by 


where ky W9 = ¥ 1-2/24,7 


WwW 
is 


conditions. 


For a liquid of specific 





2 
Ap = % W 
he? 2 


gA 


Regrouping terms 


Ap 
Fae) 


. (in*) (11) 
2 
(12) 


(in”*) 


will be defined as the equivalent flow rate referred to static 


* 
Unless otherwise stated w refers to total conditions. 


weight 7 (1b/£t>) equation (6) becomes 


(13) 


(14) 
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Determination of ky 9 by air tests of amodel makes it possible to 
predict pressure losses in a unit for particular fluids at a given flow 
rate, 

Prior to the application of this theory to the ARES model it was 
decided to construct an arbitrarily shaped flow passage which could be 
tested with air and water varying Re and M. This passage, which is 
shown in Fig. 22, has both friction and momentum losses. It was de- 
signed to operate in the same Re range for both fluids. Its inlet and 
discharge areas are equal to 4 sq. in. At the inlet is fitted a con- 
toured nozzle to minimize entrance losses into the first section of the 
pipe. Internally the model has a sudden expansion, a flow straightener, 
and a sudden contraction. The minimum flow area of 2.23 sq. in. occurs 
in the flow straightener. 

Air was obtained from the installation previously described. The 
model was attached to plenum 2, as shown in Fig. 23. Inlet total pres- 
sure is measured by the Kiel probe, and inlet temperature was taken as 
that measured in plenum 2. Discharge pressure P was measured by 
one of the probes shown in Fig. 17. Since eal tee at the discharge 
were found to be nearly uniform, except in close proximity to the walls 
of the pipe, the average of five traverse points was taken as P 


dis 


Static pressure P was taken to be atmospheric. 


dis 


Values of k, » inlet to discharge pressure loss coefficient, 


in-dis 
were established for Re (based on d. between 1.58 x 10° and 3.67 x 
10° and are plotted in Fig. 24 with the curve giving the best fit 


through the data points. Variation of k is from -9.042 to 


in-dis 
+0.057 ae Higher values of Re could not be reached due to a lim- 


iting manometer height. The highest Mach number M, obtained at the 


2] 


discharge as determined by the static to total pressure ratio, was 
0.28. A slightly higher M was undoubtably present in the flow 
Straightener. 

Figure 25 shows the test installation for the water tests of this 
model. At the 4 in. discharge of a fire hydrant wre available about 
400 gal. H,0/min. at a pressure of 50 psia. After steady flow was 
reached, the discharge was deflected to fill a 50 gallon drum. By 
weighing the drum before and after each run and knowing the elapsed 
time, the flow rate was found. Inlet static pressure Pin was measured 
just after the entrance nozzle by a 50 psig. pressure gage. 

Values of K a-dis were calculated for Re between 1.33 x 10? and 
4.88 x 10° for 4 separate runs. The method of data reduction is given 
in Appendix B, and the results are shown in Fig. 24. Except for a few 
points there is excellent agreement between water data and the curve 
established by the air test. A maximum deviation of 10 per cent is 
noted while most points fall within 2 per cent of the curve obtained 
from the air tests. At both extremes of the Re spectrum the pressure 
ies tended to fluctuate, which accounts for some of the scatter in 
these areas. What occurs at higher Re is hard to predict although it 


appears that k is tending to approach some constant value. 


in-dis 
Because of limitations imposed by the method used to collect and weigh 
water, the higher Re range could not be tested. If a flow measuring 
orifice were installed upstream of the model inlet so that high flow 
rates could be measured accurately, data at high Re could be obtained. 
These tests prove that the use of equations (12) and (15) for the 


evaluation of the pressure loss coefficient for liquid flows from air 


tests is justified. 
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SECTION 5 


DESCRIPTION OF TESTS 


In addition to the two configurations of the ARES model and the 
flow channel mentioned previously, tests were made on a full scale 
model of a nozzle tube to determine a representative loss coefficient. 
Moreover, separate tests were made with a geometrically similar model 
of a tube section in the ARES model from plane W to plane X, simu- 
lating the nozzle tubes. These tests are described in the sequence 
they were made. Table I gives an index of the runs conducted. All raw 
data from these tests are on file with Professor Vavra at the Pro- 
pulsion Laboratory. 

5.1 Oxidizer Return Sub-assembly 

In testing the oxidizer return sub-assembly, a total of six 
runs was carried out at Re between 0.868 x 10? and 1.343 x 10° referred 
to conditions just upstream of the flow measuring orifice. The maximum 
Re attainable was restricted by the height (250 cm) of manometer bank 1. 
As a result of Vavra's preliminary tests it was determined that a new 
discharge total probe assembly should be designed. Runs R1 and R2 
were conducted without taking any measurements in discharge plane I. 

The new probe assembly shown in Fig. 17 was built, calibrated, and 

used first during run R3 where alignment and ease of adjustment were 
investigated. More complete data in plane I were taken during runs 

R4, R5 and R6. At each selected peripheral location the probes were 
aligned individually with the injectors. Although the readings were 
steady and accurate, the procedure was time consuming. It was therefore 


decided to use data from every third injector tube only. 
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Pressure taps were made in the annulus by drilling a 1/16 in. 
hole partially tapped with a 10-32 thread. If the portable probe show 
in Fig. 19 was not inserted, these holes were plugged by screws. 
Because of wear on the threads in the plastic, conditions in the annu- 
lus were determined only for runs R3 through R6. 

Pressures from plenum 2 downstream through plane H were measured 
on manometer bank 1 referenced to atmosphere. Pressure in plane I was 
measured on bank 3,also referenced to atmosphere. Because of the 
large number of tubes and the slight pressure fluctuations noticed dur- 
ing testing, Polaroid pictures were taken of bank 1 for more accurate 
data recording. 

At locations 5, 6, 7 and 8 in plane H the pressures were observed 
to fluctuate considerably, especially at higher flow rates. The area 
of this unsteadiness is visible in Fig. 10. For runs Rl and R2 the 
fluctuation was dampened by clamping the pressure lines leading from. 
these taps to the manometer bank. During runs R3 through R6 the clamps 
were removed and an average was taken between the maximum and minimum 
fluctuation. This method produced the same readings as with the clamps, 
but the fluctuations could now be observed. 

The general procedure for conducting each run was as follows: 

1. The micromanometers and thermocouples were Zeroed and 
atmospheric pressure recorded. 

2. The desired pressure in plenum 2 was set and the flow was 
allowed to stabilize. 

3. Flow measuring data were recorded and a photograph of bank 1 
was taken. 

4. Plane I and annulus data were recorded. 


5. Part 3 was repeated. 
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Because of the interesting velocity distribution in the discharge 
plane and the large pressure fluctuations at locations in plane H, 
attempts were made to observe the flow in the annulus visually. Dif- 
ferent methods to introduce smoke had negative results. The first 
attempt used a 50 gallon drum within which a standard Navy MK-5 day 
signal flare was ignited. This particular flare produces about 15 
minutes of a dense, slightly toxic, white smoke. The drum was con- 
nected directly to plenum 2 and pressurized so that the smoke would go 
directly into the plenum. At the discharge of the model there was 
arranged a device to pipe the flow out of the building. It was found 
that too much smoke was produced, and the mixing in plenum. 2 was so com- 
plete that there was no streamline definition in the annulus. The 
flare produced enough heat to melt the rubber pressure seals in the 
drum, allowing toxic smoke to enter the room. For these reasons this 
attempt was abandoned in favor of introducing: smoke directly into the 
annulus. Cigars soaked in mineral oil and paraffin were ignited ina 
small container which was pressurized. Lines from this container led 
to eight 1/16 in. diameter holes spaced around the entrance to the 
annulus. The smoke partially condensed in the lines and that part which 
entered the annulus diffused quite rapidly. Although streamlines were 
faintly discernible, they were of insufficient definition to observe the 
flow. 

This latter method left a large residue of tar and nicotine in the 
model. With the air passing through water was introduced into the 
annulus to clean out this residue. It became evident then that water 
would work as a visualization fluid. The set-up shown in Fig. 21 con- 


sists of a manifold connected to a plastic liquid soap bottle. Water 
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colored by food dye was forced into the annulus through the top four 
holes by squeezing the soap bottle. This method was very successful 
and produced a flow picture of photographic quality. Motion pictures 
of runs R3 through R6 were taken using a Kodak Cine-Special II, 16 mm 
camera. Both black and white and color film were used at 64 frames per 
second. Runs R3, R4& and R5 were photographed in black and white using 
TRI-EX, ASA 120 film and a red filter. KODACHROME II, ASA 25, was 


used on runs R5 and R6. 


5.2 Nozzle Tubes 
The nozzle tubes were tested in the configuration shown in 
Fig. 26. Water was supplied by a 1.0 in. hose from a faucet where the 
pressure was higher than 50 psig. A total of 22 data podnts during two 
runs was taken between Re of 1.844 x 10° and 6.874 x 10", referred to 
inlet diameter. 
Static pressure was measured at the tube inlet by a 50 psig pres- 
Sure gage accurate to within + 0.1 psi. The flow discharged to atmo- 
Sphere. For each data point the inlet pressure was set and the flow 
allowed to stabilize. Water was then collected in a five gallon bucket 
during a measured time. Data recorded were inlet static pressure, 
weight of water, and elapsed time. The maximum Re was obtained with 
the water supply fully open. 
5.3 Geometrically Similar Model of Plane W to Plane X 
A single tube from plane W to plane X was scaled up by a fac- 
tor of 3.5 and tested with water by using the set-up shown in Figs. 27 
and 28. Water was supplied through a 1 in. hose to a 24 in. length of 
6 in. diameter pipe, at the bottom of which is installed a contoured 
nozzle to ensure smooth entry into the 1.05 in. diameter pipe section. 


For runs Wl and W2, the orifices A and B had diameters of 0.761 and 
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0.757 in., respectively. Prior to run W3 new orifices, each having a 
diameter of 0.714 in., were installed. A constant water level in the 
reservoir was maintained at all times, and the flow rate was con- 
trolled by regulating the water supply and the discharge valve. The 
Reynolds number Re, referred to inlet diameter, could be changed from 
0.415 x 10° to 4.078 x 10° during the tests. The runs were conducted 
in the same manner as those for the nozzle tubes. The pressures at 
the 6 locations indicated in Fig. 27 were read from a manometer board 
referenced to atmosphere. 
5.4 Discharge and Return Assembly 

Based on the information from these water tests two new ori- 
fice plates, with holes of 0.204 in: . diameter, were installed in the 
ARES model. A total of five runs, R7 through R11, was carried out at 
Re, referred to the 4 in. supply pipe, between 0.917 x 10° and 1.32 x 
to" Pressures from plenum 2 downstream through plane B were measured 
on bank 1 referenced to plenum 2. Bank 2, referenced to atmosphere, 
measured pressures from plane C through plane H. Bank 3 was used again 
for plane I. The test procedures were similar to those used for the 
return sub-assembly. Two cameras were used to photograph manometer 
banks 1 and 2 simultaneously. No data was taken within the second 


annulus. Total pressure in plane Z was measured by the removable probe. 
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SECTION 6 


RESULTS AND DISCUSSION 


The oxidizer return sub-assembly was tested separately, mainly to 
investigate the velocity distribution in the discharge plane and the 
stability of the flow inside the return annulus. To a large extent 
these criteria determine the feasibility of the flow passage design. 
Assuming that fuel is pumped into the combustion chamber at plane I and 
distributed uniformly around the periphery, the importance of a uniform 
oxidizer distribution can be seen. The mixture of fuel and oxidizer 
will be rich or lean depending on whether the oxidizer flow is defi- 
cient or excessive. If the oxidizer discharge varies either with time 
or peripherally, and if it is not possible to meter the fuel accord- 
ingly, a non-uniform combustion will result. This condition is unde- 
Sirable for two reasons. First, hot spots will result in the combus- 
tion chamber inducing thermal stresses in the combustion liner. Second, 
and more important, is the effect on turbine inlet conditions. The 
combustion chamber length is short, and non-uniform temperatures would 
have insufficient time to adjust prior to reaching the turbine inlet. 

A turbine blade then would pass through a non-uniform temperature field. 
An off-design condition would result and could have a detrimental effect 
on turbine performance. A major portion of the regenerative cooling in 
the ARES system occurs within the return annulus. The convective heat 
transfer in this annulus is totally dependent on the oxidizer flow 
pattern. 

Since the return sub-assembly differs from the full model only 


upstream of plane E, it seems reasonable to assume that flow downstream 
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of this point will represent the conditions that will exist in the 
fully assembled model. Actually it was found that inlet conditions at 
plane E were essentially the same in both test configurations. 

Data reduction for runs Rl through R6 is contained in Appendix C, 
part 2. The equivalent flow rate vo referred to conditions at the 
flow measuring nozzle, and the values of Re just upstream of this point 
are listed in Table II. Appendix A gives the derivation of the flow 
measuring equations. Actual flow rate w varies from 0.283 1b /sec in 
run R3 to 0.449 1b _/sec in run R6. Also listed in Table II are average 
total pressures and equivalent flow rates at intermediate planes D, F 
and G found according to equation (C2). Losses between planes will be 
referred to one of these average inlet conditions. 

It is important that the flow is distributed evenly around the 
periphery as it enters the model. In this case the procedure of aver- 
aging conditions at each plane can be justified. From water tests con- 
ducted by ee: on a geometrically similar model of a tube from plane 
D to plane X, an orifice discharge coefficient K was obtained. Vavra's 
model scaled the orifice size by a factor of 3.5 from 0.216 in. to 
0.757 in. This discharge coefficient, analagous to “Xf of equation 
(A3), takes into account some of the difference between measured and 
actual flow rate. Figure 29 shows K plotted as a function of Re 
referred to tube inlet conditions. Values of K in Fig. 29 appear some- 
what higher than anes listed by ASME? for a sharp edge orifice with 
6 = 0.7209. This is quite reasonable when the different inlet con- 
ditions are considered. The ASME data are dependent on a minimum 
length of straight pipe upstream of the orifice. For Vavra's tests an 


inlet nozzle was present only ahead of the orifice. In later tests of 
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the same orifice with some length of pipe upstream, values of K were 
obtained which agreed more closely to those listed by the ASME. 
The flow at plane D was found to be stagnated, and the pressure 


Po = P was uniform. At each measuring point in plane E an equiva- 
D 


* 
lent flow rate wy was found according to equations (C5) and (C6). 
i 


* 
Assuming that Wp existed through each of the 52 orifices the peri- 
A: 
pheral flow distribution can be expressed by a factor F. > where 


(52) wy 


5 ly 
F,=— (15) 
“D 


The total equivalent flow rate “| is given in Table II. Tables III 
and LV list values of F. in addition to intermediate quantities used 

in the calculation. Figure3l is a plot of F, as a function of peri- 
pheral location for Rl through R6. Looking downstream the angle 

is measured counterclockwise with 9 = 0 deg. corresponding to the 
three o'clock position. In each instance the flow distribution is 
almost identical. The maximum deviation of this average F. from 

F. = 1.0 occurs in run R3 where the lowest Re was encountered. [It was 
observed that most of the questionable data obtained during these tests 
occurred at the lower Re where the measuring accuracy is lower. Gen- 
erally each curve has two peaks and two valleys indicating a deficiency 
of flow at © = 110 deg. and 300 deg. and an excess at © = 0 deg. and 
225 deg. This is a relatively symmetric condition with the maxima 
appearing about 180 deg. apart. Because of this somewhat regular 
occurrence the values of F, at © = 90 deg. for runs R3 through R6 are 


questionable. Not be be overlooked, however, is the fact that the 


maximum deviation from F, = 1.0 is only 4 per cent. This fact tends to 
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establish the accuracy of Fig. 29 and to substantiate that flow enters 
the oxidizer return sub-assembly uniformly. 


Figure 30, also plotted from Vavra's water tests, shows $y» 


the loss in total pressure from plane D to plane X divided by the static 
pressure change across the orifice, as a function of Re. The total pres- 


sure Pe at plane X is found according to equation (C8). With this 
X 
pressure and the inlet conditions at plane D, the coefficient kK -x was 


calculated for each peripheral location. An average value ic for 


each run was found which is shown in Tables ITI and IV. It should be 
noted that at any peripheral location the value of kK -x differs from 
ae by not more than five per cent. 

If, as was hypothesized, the losses in this section are indepen- 
dent of Re, then ky -x should not vary from point to point. This may be 
shown 1 referring to equation (12) and recognizing that both AP), 


# , = 
and ", depend essentially on AP) The distribution of Kp -x/*p-x 


follows the same pattern as that established for F 


= 


L This is shown in 
* 
Fig. 31 for Rl and R2 and holds for the other runs also. Since w. is 


D 
proportional to Re in accordance with equation (C6), the loss amis. 
cient is also proportional to Re. 

Experimentally obtained values of the friction factor proportional 
to k as a function of Re and relative roughness are shown in Fig. 32 
for two different conditions of flow in pipes. These curves were taken 
from sei, tae For the artificially roughened pipe of Fig. 32(a) 
the friction factor increases with Re to approach a constant value, 
whereas for the commercial pipe of Fig. 32(b) the factor decreases with 


Re toward a constant value. In the former case the pipes tested were 


lined with protuberances which penetrated completely through the 
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laminar sub-layer of the boundary layer. A large part of the friction 
is then due to a form drag which can be thought of as a momentum loss. 
In the second case skin-friction alone accounts for the friction fac- 
tor. As relative roughness is increased in both cases the friction 
factor reaches a constant value at a lower Re. Realizing that in the 
present model test each tube from plane D to plane X is machined to at 
least the equivalent of commercial roughness and that ky -x increases 
slightly with Re, it can be concluded that the momentum loss created by 
the orifice predominates. Tables III and IV show values of Re corres- 
ponding to the general transition range of Fig. 32. The loss coeffi- 
cient ky -x would become constant if either more momentum loss were in- 
troduced or if Re were increased. 

Incompressible flow is a basic assumption of the theory proposed 
for this analysis. It seems reasonable to assume that the highest Mach 
number M in the model exists at the minimum flow area which occurs at 
the vena-contracta downstream of the orifice plate. At plane E an 


average area restriction factor was found using equations (C9) 


E 
through (C13). Also calculated were average values of Fee Pe and M,. 
Table V lists these values for runs Rl through R6. Applying the values 
of p E to a flow area at plane E of 3.676 sq. in. gives an effective 
flow area varying from 2.081 sq. in. to 2.161 sq. in. The next 
largest area in the model is the 3.138 sq. in. area at plane I, where 
an average M, was calculated and is listed in Table V. For all runs 
the highest M in the model occurs at plane E, reaching a maximum value 


of 0.34 in run R6. Hence it is permissible to assume that the flow is 


essentially incompressible. 
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Table V also lists values of average loss coefficients from planes 


D to F, F to G, G to I and X to F. The latter loss is given by 


“eh Der Dek oo 


From Fig. 30 are obtained values of f which can be used to pre- 
dict a » and hence os whereas the values k# are calculated from 


X 
test data. The magnitude of k is small because of the small dis- 


X-F 
tance between the two planes. [It can be said, therefore, that the 
method used to establish Fig. 30 is a valid one. 

The peripheral distribution of le is shown in Table VI. 
Except for a few isolated values, Ky -F is contained within 2 per cent 
of yep Referring to the previous discussion, concerning the rela- 
tionship between equivalent flow rate and loss coefficient, it can be 
assumed that the flow is distributed uniformly at plane F. 

The situation is different in plane G. Figure 33 shows the per- 


ipheral distribution of Kc! “ec: 


was damaged, a curve was faired in to give an estimated value at this 


Since the pressure tap at O=0 deg. 


location. Except for runs R3, at low Re, a definite pattern can be 
observed. A minimum loss coefficient occurs at 90 = 180 deg. This in- 
dicates an irregularity of the flow at this peripheral location. This 
situation can be explained from the geometry of the flow passages. 
Figure 34 shows chesatereeeeien of the two cylinders with the posi- 
tions of OQ in plane G indicated. At 9 = 90 deg. and 270 deg. the 
flow has an easier path to follow, therefore one would expect a higher 
flow rate at these points and hence a higher value of loss coefficient. 
However, at 9 = 0 deg. the flow enters the return annulus on top of 


the holes in plane H. It will be shown than the primary injectors at 
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this location pass a higher percentage of the discharge flow. It is 
entirely reasonable then that the flow is deficient as it enters the 


annulus from @ = 135 deg. to © = 225 deg. 


Table VII lists values of loss coefficient Koy for each location 
in plane H and ke for runs Rl through R6. Also shown in Table VII 
are values of oh for runs R3 through R6, found by 
‘Te Came® =Gan (17) 


Shown in Fig. 35 are values of average loss coefficients between 
different stations of the model as a function of Re. The increase with 
Re of a has already been explained, however, the maximum and minimum 
values are only 92 and 104 per cent of the average value, respectively. 
The highest value which as approaches at increased Re sts probably 
around 0.125. If loss coefficients in the low Reynolds number range 
are not .considered,.the average losses from plane X throughout the re- 
mainder of the model are remarkably constant. 

Velocity distribution in the discharge plane is found according 
to equations (C15) and (C16). Figure 36 shows a drawing of plane I 
from the discharge side and the relative position of the pressure taps 
in plane H. Table VIII lists values of v,/V for each of the surveyed 
locations for runs R3 through R6. An average of positions A, B and C 
was calculated at each location lthrough 57. This average velocity 
distribution is plotted in Figs. 37 through 40. The velocity distribu- 
tion and hence the flow distribution is essentially independent of Re. 
The excessive flow from locations 54 through 10 was anticipated from 
the conditions existing at OQ = 0 deg. in plane G. The erratic con- 


dition from locations 24 through 36 will be explained later. It is not 
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understood why the flow is deficient around location 48. Data at the 
seven locations within the annulus obtained from the removable probe 
were entirely unreliable. The model had been dropped causing damage to 
the plastic threads and making it virtually impossible to prevent 
leakage. 

The flow within the return annulus can best be described by re- 
ferring to the results of the flow visualization efforts previously de- 
scribed. It has been shown that the flow in plane G is symmetrical 
about the 6 = 0 deg. and O = 180 deg. axis. Therefore, the annulus 
is naturally split into a top and bottom half in each of which the same 
flow pattern should exist. Water introduced in the lower half had a 
tendency to pool at the lowest point because of gravity effects, whereas 
all the water introduced in the top half was carried out of the model. 
For these reasons and because of problems associated with photographing 
the flow conditions, it was decided to investigate only the upper half 
of the annulus. Figure 41 is a drawing of that half annulus, developed 
into a plane, showing the locations where water was introduced. In 
addition, data obtained by vara” are shown that give flow directions 
at locations within the annulus. These flow directions are represented 
by arrows in Fig. 41. Figure 42(a) through (h) and Fig. 43(a) through 
(h) are photographs taken of runs R4 and R5, respectively. The upper 
half of the annulus is shown as the upper 3/4 of each picture. These 
pictures were memes single frames of 16mm movie film taken at 64 
frames per second. Each sequence spans one second of operation with a 
time between each pattern of 0.125 sec. Figure 42 is typical of the 
flow pattern at Re = 1.028 x 10° and lower whereas Fig. 43 is repre- 


5 
sentative of the conditions at Re = 1.177 x 10 and higher. Several 
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similarities exist at all Re. Most apparent is the extreme instability 
in the region of plane H where the flows from the top and bottom halves 
of the annulus mix. This region is exactly that in which the pressure 
taps H5, H6, H7 and H8 are located. 

It should be noted that this is the only region in the model where 
the pressures were unsteady. Several vortices are present in this gen- 
eral area. The pattern is apparently being triggered by a huge separa- 
tion from the turning vane near H5. Although not entirely obvious from 
the sequence of pictures, it was observed that these vortices went 
through a periodic variation in both intensity and location. At best, 
it can be said that the flow in this area is thoroughly confused and 
that this situation is most likely the cause of the cyclic velocity dis- 
tribution in plane I which is evident from Figs. 37 through 40. It has 
been noted also that a deficiency of flow exists at the far right of 
the annulus. This portion of the annulus is supplied by the 9 = 180 
deg. location of plane G where the deficiency of flow has already been 
observed. The most striking difference between Fig. 42 and Fig. 43 is 
the appearance in the latter of a huge vortex in the upper right hand 
portion of the annulus. It is believed that this vortex is caused by a 
flow separation from the turning vane immediately upstream of it. A 
100 ft. roll of color film taken of runs R5 and R6 was left intact, and 
is on file with Professor Vavra. 

From the above discussion and the flow visualization studies, it 
is possible to make a qualitative design evaluation of this portion of 
the ARES assembly. The test data at the higher Re will best simulate 
actual operation, although this distinction appears necessary only when 


considering conditions within the annulus. Heat transfer will most 
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certainly be hampered where vortices exist, since the same flow parti- 
cles will tend to circulate in one spot instead of passing through the 
model. Hot spots could result in these regions. The deficiency of 
flow in the region farthest from plane H would have the same effect on 
heat transfer, only to a lesser extent. Elimination of these flow in- 
stabilities, in addition to improving heat transfer characteristics, 
would serve to reduce the losses from plane G to H and to distribute 
the discharge velocity more uniformly. One method of accomplishing 
this would be to redesign the turning vanes, making their shape conform 
more to airfoil cross sections cambered in the direction in which the 
flow is to be turned. For stress considerations, the contact surface 
area with the outer annulus would have to be maintained. Another 
scheme would be to replace the relatively large number of turning vanes 
with a few long flow dividers. It is felt that a large powtion of the 
instability occurs as a result of the mixing of flows from each half of 
the annulus, and for this reason the arranging of a flow splitter is 
recommended. 

The loss coefficients calculated during the tests of the return 
sub-assembly were duplicated when the full model was tested and will be 
discussed later in this report. Before proceeding with the tests of 
the entire unit, it was decided to determine how closely the pressure 
drop in the model, from plane W to plane X, simulated that in the 
actual nozzle aie Results of tests on the nozzle tube are listed in 
Table IX. The method of data reduction is shown in Appendix C, part 3. 


Loss coefficient 5 is plotted as a function of Re in Fig. 44. 


in-dis 
Fluctuation of the pressure gage was as much as + 0.5 psig. At low 


values of Re, where pressures were of the order of 3.5 psig., this 
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fluctuation could cause the computed % to vary between 11.1 and 8.5. 
Therefore, the dashed portion of the curve of Fig. 44 is questionable. 
The decrease of 5 with Re is quite evident, and by applying the 
rationale previously described, it seems that skin friction accounts 
for the major portion of the loss. Because of the scale chosen for t 
in Fig. 44 its variation appears magnified. Actually g only deviates 
+ 6 per cent from an average value of 10.3. The range of Re encountered 
falls well within the transition range for pipes as shown in Fig. 32. 
Data from tests on the geometrically similar section of plane W to 
plane X were reduced by the method outlined in Appendix C, part 4. 
These tests were conducted similarly to those run by Vavra for which the 
validity has been established. Table X lists values of fe rR 1-6 
and the orifice discharge coefficients, in addition to other quantities 
of interest. Points 1 and 6 refer to planes W and X respectively. 
Figure 45 shows f 1-6 and eae for runs Wl and W2, plotted as a func- 
tion of Re. The loss coefficient - differs. by almost a factor of 
2 from . ‘tere of the nozzle tubes. This condition is unsatisfactory 
because the largest loss in the entire assembly occurs in these tubes 
and the water model with orifice plates of d = 0.757 in. and 0.761 in. 
gave half of this loss only. It is also interesting to compare the 
values of K (0.757) with those obtained by Vavra. The values differed 
on an average by almost 6 to 7 per cent. The difference in the dis- 
charge coefficient is due to the influence of the inlet conditions. 
Table XI also lists the results of run W3 with orifice diameters 
of 0.714 in., whereas Fig. 46 shows E 1-6 and fener as functions of 
Re. It appears as if this orifice diameter more nearly simulates con- 


ditions in the nozzle tube. From the test of the fully assembled model, 
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values of & w-x were computed as outlined in Appendix C, part 5, 
equations (C32) to (C36). Figure 47 shows these values together with 
the curve obtained from the water tests of the nozzle tubes. The dif- 
ference in inlet conditions to plane W when the section is installed in 
the model could account for the differences between . W-xX and ic 1-6" 
Apparently though the orifice size of 0.714 in. does the best job of 
simulation. 

Tests of the complete assembly of Fig. 4 were conducted mainly to 
establish values of loss coefficients without making pressure or veloc- 
ity surveys. The methods of data reduction are similar to those used 
for the return sub-assembly and are contained in Appendix C, part 5. 
Table XII gives the data for the flow rate determined, together with 
average total pressures and equivalent flow rates in planes Z, A, F 
and G. Data at plane Y indicated a lower total pressure than that 
occurring in plane A. This impossible situation must be due to inac- 
curate measurement in plane Y. These data were therefore disregarded. 
Average loss coefficients from planes Z to A, A to F, F to G, G to H, 
H to I and W to X are listed in Table XII for runs R7 through Rll. 

The range of Re encountered in the 4 inch supply pipe, and in the 
model at plane E, was the same as for the return sub-assembly alone as 
was the value of discharge Mach number M Arranging the discharge 
sub-assembly ahead of the return sub-assembly does not disturb the 
peripheral Select ty ai see neuen that was observed in the return sub- 


assembly alone. Table VI lists values of ky : p/k oy which differ 
from a value of 1.0 by a maximum of 0.005. The distribution of loss 


coefficients from plane F to plane G is shown in Fig. 48. It is almost 


identical to that from runs Rl through R6. Velocity distribution in 
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Plane I was also found to agree precisely with that determined pre- 
viously. From the foregoing it is reasonable to conclude that the dis- 
cussion of results for the return sub-assembly is applicable to the 
entire model. 

Figure 49 shows the average loss coefficients as a function of 
Re. The scatter of eee Occurs because the total conditions at plane 
A were measured at only one location, hence a true average could not be 
obtained. However, this condition is not serious because this loss is 
only a small percentage of the total loss. The total loss coefficient 
es is obtained by adding the coefficients between intermediate planes. 
If only the conditions in planes Z and I are considered, a value io 
is obtained which differs from ee, - As Re is increased ee 
goes from 0.974 to 0.948. Evidently there exists a slight effect of 
compressibility. The difference in the two loss coefficients may be 
shown mathematically. Assume, for example, that the loss coefficient 
between planes Z and F is to be determined. The overall coefficient 


may be expressed as 


, 
— 1 ee =n ome ox 
k Fn Ce ere Er.) (18) 
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Since total pressure must decrease from plane to plane, the pressure 
ratio P. /P_ must be larger than unity. Hence, equation (20) shows 

Z A 
that Ky oF will always be larger than ky eg The total loss coefficient 


from plane Z to plane I can be expressed by 


ee 
7 alien AT eagle tll Tales ser a) 
where 
=e — 
K AoE = = SAF (etc.) 
“A 


Equation (21) was used to compute K, oe fer runs R6 through Rll. The 


values were found to be identical to those obtained from 


’ oP, 
mee ET — *2 
1 = (we ) (22) 
- 
Z 


Consider the section from plane G to plane H only. Figures 35 
and 49 show that the loss coefficient as a function of Re, if based 
only on inlet conditions at plane G, is the same for both the return 
Sub-assembly and the fully assembled model. If the method of equation 


<= ) : 
(21) were applied, it can be seen that k would increase as more 


G-H 
sections were added to the model upstream, because for the same inlet 
* 
conditions at plane G, namely, P andw , the pressure ratio P. /P 
we G tc Ca 


would increase. This method actually refers all loss coefficients 
back to the initial inlet conditions and assumes incompressible flow 
from, say, plane Z to plane H. Using re to represent the loss, 


effectively cancels any errors introduced by the incompressible assump- 


tion for flows upstream. Only in that portioncaf the model from 
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G to H is incompressible flow assumed. For this reason, it was decided 
to use the above-mentioned simple method for computing total loss 
coefficients. As an extension of this reasoning, the writer believes 
that more accurate estimates of loss coefficient could be obtained by 
adding more intermediate measuring planes. It is questionable whether 
this complication is justified since the maximum difference between 
i and — is only 5.2 per cent. 

From Fig. 49 the coefficient as appears to approach a value of 
about 0.62 of which 66 per cent is from plane W to plane X. As an 
example it will be assumed that the flow rate of the oxidizer w is 
100 1b _/sec and that its specific weight Le is 75 lb/cu. ft. The 
total pressure drop re in the oxidizer passage from the pump dis- 
charge to the entrance into the primary combustion chamber is then from 


equation (14) 


2 
_ w (144 eo: 
AP a 7 a (1b/in’) (23) 


where: ee = 0.62 ine 


g = 32.17 ae Seca 


Or 


2 
> = 0.62 (100)" (144) 2 379 aByin? (24) 


“2-1  (@2m07) =075) 
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SECTION 7 


CONCLUSIONS 


The average loss coefficient from a position just downstream of 
the oxidizer pump to the discharge of the primary injectors is 0.62. 
This experimentally determined coefficient should remain nearly con- 
Stant near the specified operating range since the anticipated Reynolds 
number is sufficiently high. The model section from plane W to plane X 
with 2 orifice plates, having holes of d = 0.204 in., represents the 
actual nozzle tubes. The pressure loss coefficient found in the air 
test of this model section agreed closely with the coefficient deter- 
mined from the water test of the tubes. The peripheral distribution 
of flow within the model is very nearly uniform until it enters the 
annulus. Here the flow becomes deficient in a 45 deg. sector centered 
on 9 = 180 deg. in plane G. This causes a deficiency of flow in that 
portion of the return annulus farthest from plane I. 

Because of compressibility effects a difference exists between 
loss coefficients determined from inlet and discharge conditions only, 
and those obtained by a summation of the intermediate loss coefficients 
between adjacent planes. The flow channel used to check the basic 
theory by water and air tests was neither designed nor instrumented with 
this difference in mind. [It would be of interest to construct a more 
complicated system of flow passages for this purpose and to investigate 
the effect of compressibility more fully. It is believed by the writer 
that more intermediate planes where exact conditions can be established 
will lead to more accurate estimations of loss coefficient. 

Due to the fact that the major portion of the pressure drop is 


caused by momentum losses, it would be inconsistent to apply standard 
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theories that relate heat transfer to pressure drop with the objective 
of evaluating the heat transfer coefficients at different locations in 
the flow passages. However, the vortices within the return annulus, as 
evidenced by the photographs of the existing flow pattern, will cer- 
tainly be detrimental to heat transfer in this region. These vortices 
could be largely eliminated by changing the design of the turning vanes 
and installing a flow divider where mixing now occurs between flows 
from the top and bottom of the annulus. 

Model tests of this type provide a natural bridge between the 
completion of a design and the building of a prototype. It is estima- 
ted that a cost factor of 20 to 1 exists between the actual hardware of 
the ARES unit and its corresponding flow model. Prediction of the 
pressure loss coefficients of complicated passages analytically is, at 
present, an impossible task. Except for the slight effects of Reynolds 
number, the experimental method of predicting pressure loss coefficients 


presented in this paper is both valid and accurate. 
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RUN NUMBER 


RMW1, RMW2, 
RMW3, RMW4 


RI, R2, R3, 
R4, R5, R6 


W1, W2 


W3 


WV1, WV2, WV3 


R7, R8, R9, R1O, R11 


TABLE I 


INDEX OF RUNS 


TEST DESCRIPTION 
Air test of flow device 


Water test of flow device 
Air test of oxidizer return sub-assembly 


Water test of geometrically similar model 
plane W to plane X orifice plates of 
d = 0.761 in. and d = 0.757 in. 


Same as Wl and W2 with orifice plates both 
of d = 0.714 in. 


Water tests of geometrically similar model 
plane W to plane X 


Air tests of full ARES assembly 


d1 










RUN 
AVERAGE OF ALL 
DATA POINTS 


OXIDIZER RETURN SUB-ASSEMBLY, INLET ORIFICE 
DATA AND EQUIVALENT FLOW RATES 


TABLE II 






; (1b/in?) | 14.768 | 14.768 | 14.778 | 14.778 | 14.792 | 14.734 
u 


P tm 
2 
P (1lb/in ) 
Ah, (in H,0) 
Oo 
di 
Oo 
ie lle 
Re(10">) 
d = 4.07" 
te 2 
yh oh 
P (1b/in*) 
D 
- 2 
P (ib/in*) 
E 
F 
ze 2 
P, (1b/in*) 
G 
* 2 
et 
voy 
% 2 


cafe fe fw |e |e 
14. 


<2 fname en 
152558 75132 
-54291 .81237 
54708 92829 


~71732 | .70137 


76914 | .74870 


a2 


TABLE III 
FLOW DISTRIBUTION THROUGH ORIFICE PLATE, 
LOSS COEFFICIENT PLANE D TO PLANE X 


Pan. * 
* 


W = .71732, K=. 


RUNS R1, R2, R3 


ae ees 2 |. 
7 








R1 " TOE. GtEp - 20D B= 1161.4 
Poe 163.10 | 159.40 | 155.05 | 151.80 | 159.50 | 163.10 | 159.10 | 158.45 
(Cm HO) 
3.503 3.468 3.426 3.391 3.468 3.503 3.466 3.458 
d=0. 
x 
YD an7) -01417 | .01403 | .01386 | .01371 | .01403 | .01417 | .01402 | .01399 


pe fice 1.016 |1.005 10.994 11.005 {1.027 |1.016 |1.014 
Reeewareiita7s |mees | ize) | aidles | 21245 | 273 e42” ies 
ave. | .1239 
ke = 
R2 wy = -70137, K = .79, Pe gues > 30 AER = 1152713 
Py Fe 151.25 | 147.80 | 143.85 | 141.70 | 147.65 | 151.60 | 147.60 | 147.33 
(Cm H,0) 
Re(io-*) 13.363 | 3.328 13.286 | 3.264 | 3.326 | 3.366 |3.326 |3.324 
d=03 
%* 
wp (in?) .01377 | .01363 | .01346 | .01337 | .01362 | .01378 | .01362 | .01361 
i 
oF, fieoa 1.010 |0.998 10.991 |1.009 |1.022 {1.009 11.009 
(in-*) |.1250 |.1221 |.1188 |.1171 |.1220 |.1252 |.1219 |.1217 
-X 
ave. | .1217 
* 
W. 


R3 = 525585 KP .7957P SP - S05 A Pa = 1095.29 


D D D-E 
Py Pe 69.70 
(Cm H,0) 








J 













X 
65.50 |68.95 |64.75 |66.85 |71.45 |67.55 |67.35 
Re (10"*) 12.249 12.187 |2.238 |2.173 |2.207 |2.276 [2.218 |2.214 
d=0.3 
* 
-01004 | .00976 | .00999 | .00970 | .00985 | .01016 | .00990 | .00988 
F, 0.994 |0.965 |0.989 |0.960 |0.975 {1.006 {0.979 |0.978 


-1128 . 1060 .1116 . 1048 . 1082 ~1156 .1093 . 1090 
avg. |.1096 
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TABLE IV 
FLOW DISTRIBUTION THROUGH ORIFICE PLATE, 
LOSS COEFFICIENT PLANE D TO PLANE X 


RUNS R4, R5, R6 





* . — 
me =) .61799, K= .79, P, =P, - 3508 O Po, = 1122.90 
104. 


mM 
=(). 


lin 11185 1.1127 par60  |2094 1150 |.1208 |.1149 {.1138 
Ky -x 
av eb 51 


0 
Ww .68534, K = .79, P. =P. - .50 A P 





= .79, P =P = 1181.06 


D 5D D-E 
PD Pe 182.90 |176.05 |176.90 {171.35 1180.90 |185.60 |177.90 1177.05 
(Cm H,0) 


. . z yg 7 1152.88 
X 
Py Pe 140.95 1135.05 |137.25 {132.05 1137.95 1142.95 1136.45 1135.95 
(Cm HP) 
Gg =U 


W. 





Re (10°*) =|3.681 |3.618 {3.626 |3.574 13.661 {3.706 {3.606 1|3.628 
A= 5 
Mt (an?) | 02479 .01454 |.01457 |.01436 |.01471 |.01489 order |.014se 
ae 1.006 h.008 |o.994 por ft.030 fi.o14 1.009 
D-X 
ave. |{.1245 
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TABLE V 





AREA RESTRICTION PLANE E, MACH # PLANE 
LOSS COEFFICIENTS PLANES D TO F, X TO F, 


os 2 
1125.75 


(cm H,0) 1169.33 1159.50. 50 1155.63 | 1184.63 
See 

mom fm fm [oe 
Fie fe fn fom fe [es 
[im | se | see | on | ow 
A a ae 
; 1161.2 1148.11 {1174.90 
lb 
ae 0293 0288 0295 
fos fn fo 


0501 |.0508 |.0517 











+ 





aS 


rg] 
rt 

jl 
7-~ 

Q 

9 

So 
“Oo 
ee” 


G-I 


oa 
N 


TABLE W 


PERIPHERAL DISTRIBUTION OF LOSS COEFFICIENT 
PLANES A TO F AND PLANES D TO F 


| 307 | c2e om || cor 
re [es oer om [1 [toe 
as [os oer om we [over 
"oe [e[ or oe | [oe 
os [a7 [sr ee) oe [ as 
[rr [or [ese] ose) ios) om 
re [oss coe] ome ose 
oe [i [or oe ee es 
ao [ae [or [oe 00 | es 
an [seo ono cnr ee 
Fa 
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DATA 
POINT 


rf 


aa 
a 
Pu 
© 
co 
™ 
+ 
a 
Hl 


atm 


P 


13 


14 





~ | * 0.693 19.194 | 27.950 | 29.411 





TABLE IX 
RESULTS OF WATER TESTS WITH SCALE 
MODEL NOZZLE TUBE 





n°/ sec. 






Ib/sec. 


0.656 18.169 | 26.000 | 27.304 


5.675 | 9.119 
30.225 31.856 
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Oo 
= 









10.728 
10.815 


TABLE X 
RESULTS OF WATER TESTS WITH GEOMETRICALLY SIMILAR 
MODEL OF SECTION BETWEEN PLANES W AND X 


W eres re K(761) | K(714) |K(.714) 
ib/sec | in3/sec | d=1.05 | !-6 1-6 |pos.B | pos.A | pos.A |pos. B 
3 












UN 


| | 


oo 


2557 le 29 5.419 


1.124 | 31.149 


76 
rae | | 
ae | 
as] 
imal i 
rus || 
ne eee al 
rss | 
cas | 
36) a | 
C40 | 
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4 


ke 270 
Desire 


rai.1es [2 


0.878 
: 9.617 


2720 
. 768 
ES 
Od 
. 780 
774 
. 780 
. 766 
767 
»/22 
743 
745 
- 746 
. 762 
. 760 
oe 
745 
754 
747 
742 


W2 


[2.357 


Ta094 [at.4o7 | 3.265 | 1.452 | 9.20 | | | 760 | 7a 
T1509 | as.020 | 3.073 | 1.007 | 9.0% | |_| 767 | 704 
1-090 | as.020 | 3.093 [ 1.0 | v.20) |_| .756 | .73_ 
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TABLE XI 
OXIDIZER DISCHARGE & RETURN ASSEMBLY, INLET 








ORIFICE DATA AND EQUIVALENT FLOW RATES 
2 
P * 14.743 1b/in 
atm 


RUN 
P (lin 18.055 | 18.943 | 19.986} 21.179] 22.220 


= 19.986 
Ah, (in H,0) 0.991 | 1.198 | 1.448 | 1.736 | 1.951 
plen 
= Ty 
Re(10™) gy gy 0.917 | 0.985 | 1.117 | 1.276 | 1.320 
2 
P, (1b/in”) 17.769 | 18.651 | 19.602| 20.750] 21.721 
Z 
P. Gnas) 17.597 | 18.388 | 19.296] 20.393| 21.330 
A 
P, (1b/in”) 15.589 | 15.827 | 16.114] 16.462] 16.752 
F 
P. any ane) 15.464] 15.668 | 15.913] 16.216] 16.463 
G 
wy ena) 50916 | .54449 | .58750| .62596| .64702 
a 9 
- 2 
wy Ge .58037 | .64164 | .71467| .78901| .83895 
Wo ca 58505 | .64817 | .72367| .80098| .85364 
P 18.005 


18.890 | 19.894; 21.093} 22.132 





— : 2 
. (1lb/in’) 
O 
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TABLE XII 
LOSS COEFFICIENTS PLANES Z TO A, A TOF, 
F TO G, G TO H, H TO I, W TO X 


DISCHARGE MACH NUMBER 








P 407.03 405.83 
22 





= f=f=t=[> 
ie (in”*) 0486 0493 

P, (1b/in’) 15.845 16.801 

-/\F 1.458 1.445 

3908 4193 


9950 10.560 | 10.698] 10.973] 11.329 









Combustion chamber 
primary 


Injectors 


NOZZLE 


FIG. 1 
ARTISTS CONCEPTION OF ae 
ARES PROPULSION SYSTEM. ie 
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FIG. 42 (con't) FLOW IN RETURN ANNULUS RUN R4 





FIG. 42 (con't) FLOW IN RETURN ANNULUS RUN R4 





FIG. 42 (con't) FLOW IN RETURN ANNULUS RUN R4 





FIG. 43 FLOW IN RETURN ANNULUS RUN R5 





FIG. 43 (com t) FLOW IN RETURN ANNULUS RUN R95 





FIG. 43 (con't) FLOW IN RETURN ANNULUS RUN R5 





FIG. 43 (con't) FLOW IN RETURN ANNULUS RUN R5 
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APPENDIX A 


DERIVATION OF FLOW MEASURING EQUATION 


The pressure drop PY = Py across the flow nozzle is referred to 


the absolute pressure Py upstream and is to enter calculations in the 


Pu-Pd 
Pu 





dimensionless form Thus any consistent unit of pressure may be 
used. 


The equation for flow rate w (1b_/sec) " is 


WwW = Cae ate 2 9 eu (Pu Fag) (Al) 
where; 
eo - nozzle discharge coefficient (independent of Re) 


Y) - expansion factor to account for compressibility effects 


Ay - area based on nozzle diameter ‘én = 7.17263 eee 

f - area multiplier to account for the thermal expansion of the 
nozzle 

f - Reynolds number correction factor 


g - 32.17 eayeee 
@€u - upstream density (1b_ /£t”) 
Pu - pressure upstream Guin.) 
2 


Pd - pressure downstream (1b/in. ) 


Using the equation of state for a perfect gas 


. * (144) 
Cu a 
R ui 


(A2) 


I 
VDI Durchfluss Messregeln. DIN 1952, VDI Verlag, Berlin, 
1943, ope. 
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where R = 53.35 ft-1b/1b °R, and T is the upstream temperature (°R), 


equation (Al) can be rearranged to give the equivalent flow rate 


Ww = wea ES. ia ie Ag iby he “ 
a VF 


The dimensionless coefficients must be investigated for the noz- 


zle with diameter ratio 6 = 0.7434 and m = B “ = 0.552 where e = d/D. 


(int) (A3) 





The discharge aEaameionta a for m * 0.552 is found to be 1.114. 


The density coefficient Yy is obtained from the ASME” standards. 








For air with ¢ = cfc, = 1.4 and 6 = .743 using x = Buske there is 
0.286 
/— a 
Y ={/ 2.433 {i= x) eles 226.) 7-428 oe) 
x /- 0.305 (1-x) 


For stainless steel the area multiplier f. can be expressed Bye 


f, = /,0 + 0.00/5 | Sie oo (A5) 


The design of the nozzle is such that r. = 1 for all pressure ratios. 


Combining the above quantities equation (A3) becomes 
x . T, 
W = Wane WK = /i=5 Y 
uUu Pp y $ 
U 


eThid. Appendix 6. 


(A6) 





Soupp lanier to ASME Power Test Codes, Chap. 4, Part 5, ASME, 
New York, 1959, p. 76. 


Stearnes, R. F. et. al. Flow Measurements with Orifice Meters, 
Van Nostrand Co., Inc., New York, 1951, p. 257. 
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APPENDIX B 


DATA REDUCTION FOR AIR AND WATER TESTS WITH FLOW CHANNEL 


B. l Air Tests 
A total of seven data points was observed. The quantities mea- 
sured were: 
Hi - pressure upstream of flow measuring orifice against atmo- 
sphere (in-Hg) 
4h, - change in static pressure across orifice (in.H,0) 


qe Oo 
ti temperature upstream of orifice ( F) 


Caen - temperature in plenum 2 (°F) 

Hage - manometer reference (cm H,0) 

HS - inlet total pressure at Kiel probe against atmosphere 
(cm H,0) 

Hea discharge total pressure against atmosphere (cm H,0) 


i - 30.01 in. Hg 


3 
The equivalent flow rate Ww, is given by equation (A6). It is 
convenient to use the model inlet as a reference condition. Then the 


equivalent flow rate is 








W = WwW (in*) (Bl) 
in uw P 
re 
where P = P. = 0.49112 } H, + 30.01} Cipyanias) 
u 
; : ; 2 
Peale Pe 0.49112 | 0.02894 (Ho - Hy) + 30.01 | Gey ini 


Reynolds number Re based on inlet diameter d is 


= Vd 4 ow 
Re a ne = ras (B2) 
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Equation (B2) can be rewritten as 


wa fate 4 ( 12) 
Ee _- | Ry 


Air temperature for these tests varied between 85 and 95°F. Introducing 





(B3) 


a representative viscosity of 3.9 x 1077 ipeeteey nea in addition to 


other non-varying quantities, into equation (B3) gives 


x 
Re = Win Fin (4.728 x 10°) (B4) 


fle, 


The loss coefficient is determined by 





P _~P 
in ta, 2 _ 74 
Kin dis " “ts fy, 8) (in a 
Pin _ 


where P = 0.49112 | 0.02894 (Hae es) + 30.01} anys) 
t atm ES. 
dis dis 


B.2 Water Tests 

A total of four runs was made, each consisting of several data 
points. The runs were numbered RMW 1, etc. The quantities measured 
were; 

Pin _ inlet static pressure Guyane gage) 

We - weight of water (lb) 

t - elapsed time (sec) 
Volume flow rate is obtained by 


\ee (100) CUne aes | (B6) 
3.609% 


Re based on inlet diameter is 
ee cee 
B/ 
Ay 774 4AV (B7) 


P21 


where = kinematic viscosity, for water at 72°F = 152.257 ial eee 
Equation (B7) is reduced to 
Re = cf 239.1645 (B8) 
Because the inlet and discharge areas are equal 


AP ™ OP a-dis a Pin 
in-dis 6 


The loss coefficient is determined by 





ie Vs /3.940% FP. 
C - x = a 
oe i, . (in) (9) 
w * w ~ 


where ¥ L specific weight of water at 72k. 
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APPENDIX C 


DATA REDUCTION FOR TESTS WITH ARES MODEL 


C.l1 General 

Methods of data reduction used to obtain loss coefficients and 
other meaningful quantities from tests on the ARES model (as described 
in Section 5 of the main body), are discussed in the following Sections. 
They are listed in the same sequence as used for the performance of the 
tests. 
C.2 Oxidizer Return Sub-assembly 

The quantities necessary to obtain mass flow rate are identical to 
those described in Appendix B. The equivalent flow rate w is given by 


equation (A6). In addition the following quantities were recorded 


ieeen - pressure in plenum 2 (cm H,0) 
H tm i atmospheric reference, manometer bank 1 (cm H,0) 
He - pressure at Kiel probe (cm H,0) 


Hy - pressure in plane D (cm H,0) 
bon - static pressure at selected locations in annulus (in .H,0) 
Hon - total pressure at selected locations in annulus (in -H,0) 


H - atmospheric reference, manometer bank 3 (in.H,0) 


atm 3 
hy - static pressure in plane I (in. H,0) 
H, - total pressure in plane I (in.H,0) 
oO 

Coton temperature in plenum 2 ( F) 
From the polaroid pictures the following quantities were read 

h, - static pressure in plane E at holes 7, 13, 20, 26, 33, 39, 
46, 52 (cm H,0) 


hy - static pressure in plane F at holes 1, 9, 17, 25 (cm H,0) 
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Ho - total pressure in plane G between holes 32-1, 4-5, 9-10, 
13-14, 17-18, 20-21, 25-26, 29-3 (cm H,0) 
hy - static pressure in plane H above discharge holes 11, 14, 18, 


22, 26, 30, 31, 34, 38, 42,°46, 50 (Com H,0) 


Re i we, a + (is) aa *% ln 0.09062 (cl) 


=- W SS 


VT WlaonfR, it 


The loss coefficient k for flow between each plane is referred to 
an inlet condition. The equivalent flow rate must then be referred to 
conditions at intermediate planes. For example,an average equivalent 


flow rate at plane F is 


er eye emits (C2) 
e aaa —_—- 
Pee Ty 


z a 2 
oe? 0.49112 § 0.02894 (Ay He) +P. (ib/in ) 


Tp - temperature in plane D assumed to be the same as that in 
plenum 2 
re and T - determined as in equation (Bl). 


—*k:-—* 
In a similar fashion the equivalent flow rates Wy > Wo are determined. 


A geometrically similar model of a tube section from plane D to 
plane X was tested by Vavra. Data from these tests plotted in Figs. 29 
and 30 include the orifice discharge coefficient K, and the loss in 
total pressure . D-x from plane D to plane X divided by the static 
pressure change across the orifice. For each orifice preceding a mea- 


suring tap in plane E, an equivalent flow rate is found by 
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o, (C3) 
L 
where Aa =0.0363 elite of orifice opening). 
E and P, are determined as was F above. For air with 
F 
8 ao /c = 1.4 from 
Pp ov 
Y = 1.0 - G@&<36 FB Pe (C4) 
t 
Fy 
Equation (C3) becomes 
ie 
W = 0,05 20a KY ieee 
, (C5) 
rf We 


Equation (B4) is modified to find Re referred to tube inlet diameter 


* P, 
NY TT 


oD 


r 
< 


?7 
Re ( 0.323% Io ) (C6) 
This relation contains a representative value of viscosity for the tem- 
perature range encountered in the tests. An iteration is carried out 
uSing equations (C5) and (C6). A value of K is assumed and equation 


* 
(C5) is solved. Using this value of Wy > Re is calculated from equa- 
i 


tion (C6). A corrected value of K is then read from Fig. 29 and the 
process continued. Generally, one iteration is required only. 
The peripheral flow distribution at the orifice plate is found 


from 


) we 
F. = (52 a," (C7) 
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Tht total pressure P_ is found from 
X 


Rot feo Sag (RR) = 


where { is read from Fig. 30 for the Re found with equation (C6). It 


is possible also to estimate Pe by using equation (C7) with { = 0.48, 
E 
a standard value for sharp edge orifices with 6 = 0.72. 


The flow area in plane E is 3.676 sq. in. It is possible to find 


an average area restriction factor d E’ The equation of continuity for 


average conditions is 


= Ww 
V= — _ (C9) 
eae 
Also 
ve a. P 
e > / ¢ 4 P, (C10) 
Combining equations (C9) and (C10) gives 
. 2 be 
W a 
ae HER PS cou 
eA ¢ 
Using the equation of state for a perfect gas P = ee equation (C11) 
becomes 
VV; RT A’ g , Cc =) 7 P i a 
oS 


my 
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Referring equation (C12) to static conditions at plane E and solving 


for d 


2 te 
2 ae’ = 
?, = L ( ea) : Pe 
2 A. FE . (C13) 


where : 





W wo / 
= © 
Ee. Oo P ae 


An average Mach number M is determined from the pressure ratio 


Bee . Assuming no change in total temperature across the orifice, 

E 
ive .5 be % Tp? it is possible to find the temperature ratio needed to 
E 


solve equation (C13). 


With the equivalent flow rates and equation (12), average loss 


coefficients ae k 


a it _— 
F-G’ Ke Keay and k,_, are determined. For example, 


the average loss coefficient from plane F to plane G is 


a 7 P = 
eit tee /pn OM Cn) (C14) 
F-@ B, (w ™) 


At a particular location in plane I the discharge velocity is pro- 


portional to the square root of the dynamic pressure. 


vie [Hate oe H | ‘ Hato s 7 h. | = Ga. (e.tig) xcus) 


At each location the velocity relative to the average is then 


fn 
Vi) = F — S Vd (16) 
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C.3 Scale Model Nozzle Tubes 
The following quantities were recorded 
t' - elapsed time (sec) 


i weight of water collected (1b) 


Pin . inlet static pressure Give a gage) 
t, ~ water temperature (°F) 

Aj, 7 inlet area (0.07843 the) 

Pela discharge area (0.11341 ee) 

bem - barometric pressure (in Hg) 


Flow rate w was determined by dividing the weight of water by the 
elapsed time. More useful for calculations is the volume flow rate 


\F = w_ (ioe) (in / sec) (C17) 
3.60945 


Velocities at inlet and discharge can be written as 


V = MAA (C18) 


Assuming atmospheric pressure at the discharge the loss in total 
pressure through the tube is 


pererr: a 1A, ea - °/q (ve v, ) (1b /in*) (C19) 


in dis 


In all tests the water temperature varied only from 55° to 59°F. For 


t = 57°F equation (C19) becomes 


-3 
a F P. + (3765 x Jo)? 


Medic Ug (C20) 


Re referred to an inlet diameter of 0.316 inches is obtained from 


5 


equation (B7) with V7 = 160.279 x 10% needs 


4- 
Re = ¥ x 10 (C21) 
oS 
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Using the inlet as reference, the loss coefficient c is 


P 
i = = Tin-dis (24 
in- dis ee oe ee 
a 


A 
ore 4 - 0.935 x /0O Eee 
4 


in 


(C22) 


C.4 Geometrically Similar Model of Section from Plane W to Plane X 





Figure 28 shows 6 manometer tubes in a schematic manner. Actually, 
Stations 2 through 6 have 2 taps each so that hy (in -H,0) is an average 
pressure at Station 2. In addition to pressures, the mass flow rate 
and water temperature were recorded as described in sub-section C.3. 
Volume flow rate is determined as in the former tests. Reynolds number 


referred to an inlet diameter of 1.05 in. is 


Re= Vx 107 


C23 
13.218 Se 


for average temperatue conditions. 

It is assumed that taps 3 and 5 are placed such that the flow has 
the same area as the orifice opening. The error in this assumption is 
absorbed in the flow coefficient K. 


Across an orifice the equation of continuity gives 
A. Veg ns Ne Vee (C24) 


where d ~- station ahead of orifice 
or - station behind the orifice 
Expressing terms in in. of liquid head, Bernouilli's equation, ignoring 


losses, gives 


LZy 


H, = Va + Ng = fleet Vor + h, (C25) 









£4 Chg - he) 


ee oe (C26) 
/- p* ; 





If losses and other factors are included in a coefficient K, the volume 


flow rate can be written as 


es w/e = KA, V,- (C27) 


Equations (C26) and (C27) are combined to give 





2 





(1- B* ) g oo, (1- A") (C28) 


ta 
Ase ap e hor A ies. 


Three orifice sizes were tested. Table C-I lists factors for each 


orifice. 


TABLE C-I 


GEOMETRIC FACTORS FOR VARIOUS ORIFICES 


4 
DIAMETER P i : 
0.757 0.7209 0.7298 0.45007 
0.761 0.7248 0.7241 0.45484 
0.714 0.6800 0.7862 0.40039 


Using Station 6 as an example and assuming stagnation conditions 


at Station 1, the head loss HL is written as 


Him, = Aon, - Me = OAL (C29) 
24 
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Compared with other losses, the loss through the inlet orifice is 


negligible. A loss coefficient from the beginning of the 1.05 in. pipe 


to Station 6 is then 


Aoi (29) 
i % ize £ (C30) 
2 
nw, 
Introducing the volume flow rate 
= ; S Hf 
Pi Ok et ew (C31) 


Ge 


C.5 Discharge and Return Assembly 


In addition to the quantities measured in the return sub-assembly 


the following pressures were read 


h - 
O 


H - 
O 


en 


static pressure at entrance of model (cm H,0) 
total pressure at entrance of model (cm H, 0) 
pas atmospheric reference, manometer bank 2 (cm H,0) 

Static pressure in plane A, 6 equally spaced taps (cm H,0) 

total pressure in plane Z (cm H,0) 

static pressure in plane A at holes 1, 14, 27, 40 (cm H,0) 
total pressure in plane A between holes 7-8, 20-21, 33-34, 

46-47 (cm H,0) 

Static pressure in plane B at holes 7, 13, 20, 26, 33, 39, 

46, 52 en H,0) 

Static pressure in plane C at holes 7, 13, 20, 26, 33, 39, 

46, 52 (cm H,0) 


- pressure in plenum 2 referenced to atmosphere(cm Hg) 
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* 


Equivalent flow rates found as in Section C.2 include Ww We 


* 


3 


a We and Wo: These flow rates are used to calculate the loss 
== oe SSS — ly ; 
coefficients A? K.-R? ke-G? Kew? an Kot 


Using Re with the 4 in. diameter of the supply pipe as character- 
istic length, the argument, K is read from Fig. 35. Total pressure 


at plane X is then given by 


- = as oe _ «2 ° 
ey 7 on 7 {ky 7 wy (C32) 
where : - as 
a is ae 7 Kx oF 


An average Re of the flow from plane W to plane X can be written as 


p 
ye a tw (0.02364) (¢33) 
(> A 
where: 
x * 
Maw SEN, Pe ween 
- (C34) 
ad sa 
and ee 
Pos Fe 
Ce bala.) 4 as (C35) 


An iteration of equations (C33), (C34), and (C35) is carried out by 
assuming an initial value of { from Fig. 46. The process converges 
rapidly because { has a limited range of values. 


a a4 = 
With the computed values of P. and Wy , the coefficient ky -x is 


calculated. From equation (12) 


2 
2A° 


Sw se = (C36) 
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